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An upper limit on the masses of scalar superpartners in split supersymmetry is 
. found by considering cosmological constraints on long-lived gluinos. Over most of 

Qh| parameter space, the most stringent constraint comes from big bang nucleosynthesis. 

A TeV mass gluino must have a lifetime of less than 100 seconds to avoid altering the 
abundances of D and 6 Li. This sets an upper limit on the supersymmetry breaking 

! scale ms of 10 9 GeV. 
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1 Introduction 



In split supersymmetry PQI2] (see also [3]), the superpartner spectrum differs drastically from 
that of traditional weak scale supersymmetry. The fermionic superpartners are present at 
the weak scale, while the scalar superpartners have much larger masses, at a scale ms- This 
represents a new scale, which is a priori undetermined. Verification of split supersymmetry 
will require not only the detection of the new TeV-mass states at colliders, but also the 
observation of indirect effects of the full supersymmetry present above ms- 

Most probes of mj are only logarithmically sensitive to this scale J3J. For example, 
couplings whose values are determined by supersymmetry deviate via renormalization group 
flow once supersymmetry is broken. The possibility of using such effects to determine ms 
in split supersymmetry was explored in |5] . There is one observable that has power-law 
sensitivity to the supersymmetry breaking scale-the gluino lifetime. Thus, measurement 
of the gluino lifetime would allow a precise determination of ms- This motivates us to 
look at cosmological constraints on the gluino lifetime. For large scalar masses, the gluino 
can become long-lived. The possibility that the successful predictions of standard big bang 
cosmology might be altered by a long-lived gluino presents an opportunity to constrain 777.5. 

The effects of the late decaying gluinos depend on their annihilation cross section, which 
determines their relic abundance. The calculation of this cross section is complicated by the 
strong interactions of the gluino [SJ Hj . The relic abundance calculation requires an estimate 
of non-perturbative annihilation processes after the QCD phase transition. We discuss the 
computation of the relic abundance of the gluinos in Section |21 

In Section El we use the results of this relic abundance calculation to place limits on 
the gluino lifetime. A particularly strong constraint comes from Big Bang Nucleosynthesis 
(BBN). We find that for TeV mass gluinos, the lifetime is limited to be less than 100 seconds. 
Finally, we relate these lifetime bounds to bounds on the scale ms- 

2 Annihilation of Gluinos 

The relic density of gluinos prior to their decay is determined by evolving the Boltzmann 
equation. This determines the freeze-out temperature where the expansion rate of the uni- 
verse is balanced against the annihilation rate of the gluinos. There are two regimes of 
annihilation that we will consider separately. The first era is before the QCD phase transi- 
tion when there are free gluinos in a QCD plasma. The second era is after the QCD phase 
transition when the gluinos have become confined in color neutral i?-hadrons. The annihi- 
lation cross section in this second period can, in principle, be much higher than in the first, 
thus leading to a second period of annihilation. 
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2.1 Perturbative 



The perturbative annihilation of gluinos is well known [6J and at low velocities the cross 
section is given by 

where iV| ff is the effective number of flavors of quarks and includes a phase space suppression 
for the top quark. We use the full expression for the annihilation cross section jH] in our 
numerical calculations. 

Because the velocities of the gluinos are small at freeze-out, there is a Sommerfeld en- 
hancement due to t-channel exchange of multiple gluons jSJ- Each loop of gluons gives a 
na s /v enhancement but can be resummed into an overall enhancement 

^iia s /v ~ira s /v , . 
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Taking the enhancement into account, the cross section becomes 

27vr 2 a 3 , t t 27 Nf 7r 2 a 3 

v tot (9~9 - 99) = T^r- 2 + 0(v-') a**Qg - qq) = / 2 2 " + O n" 1 . (3) 

Because the Sommerfeld enhancement is a long-distance effect, one might worry that this 
effect is suppressed due to thermal effects from the QCD plasma. However, the effective mass 
of a gluon in a plasma is a s T, while the typical momentum associated with the Sommerfeld 
enhancement loops is a s rrigV ~ a s ^rrigT ^> a s T. Thus, the Sommerfeld enhancement 
remains effective. 

Sommerfeld enhancement can be interpreted as annihilation through off-shell bound 
states of gluinos. One could also ask whether bound states are formed on-shell. This requires 
radiating away energy. This will suppress the formation of such states by perturbative powers 
of a. We will revisit bound state formation in a non-perturbative regime in the next section. 



2.2 Non- Perturbative 

After the QCD phase transition, the gluinos hadronize and form R-hadrons. Very little is 
certain about the spectroscopy, quantum numbers, or couplings of these strongly interacting 
particles. In order to have a reliable calculation of the relic density of gluinos, we must 
nevertheless estimate the annihilation cross section of R-hadrons. Estimates of this cross 
section have varied widely, due to different assumptions about the relevant hadronic physics. 
A commonly made assumption is that the annihilation cross section goes as ttR 2 ~ l/w^, 
where R is the geometric size of the R-hadron. This cross section would be roughly dozens of 
millibarns, much larger than the perturbative annihilation cross section. In this section we 
argue that this is not the case. There is a huge separation of scales between the gluino mass 
and the hadronization scale. This makes it hard to exchange momentum efficiently between 
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the QCD cloud and the gluino core. Thus, we expect the light QCD degrees of freedom to 
decouple from the annihilation process. 

Geometric (a ~ R 2 ) cross sections for i?-hadron scattering are applicable in the case of 
q 2 < Aq CD . Such processes do not probe the interior structure of the hadron because they 
transfer low momentum, and the hadron appears like a solid object. On the other hand, 
gluino annihilation is a large momentum transfer process and probes the interior structure 
of the hadron. The partonic picture of the hadron becomes relevant, and we expect the 
annihilation cross section is set by the size of the parton, a ann ~ m g~ 2 - 

To explore this further, we can rely on simple quantum mechanical scattering arguments. 
Despite the strong dynamics, the inelastic .R-hadron cross section must still satisfy partial 
wave unitarity 

*^<_J_(2j + l). (4) 
J m~v 2 

Notice that the cross section scales as the de Broglie wavelength of the gluino rather than a 
geometric quantity. However, this does not prove that the annihilation cross section scales 
as (rrigv)~ 2 . If many partial wave cross sections contribute, they can add up to a geometric 
cross section. For example, if all partial waves up to j max = rrigvR contribute, where R is 
the size of the R-hadron, then Eqn. H]sums to give a geometric cross section, R 2 . For this 
to occur, high angular momentum partial waves must contribute (i.e. j ~ 30 for a TeV scale 
gluino at temperatures around Aqcd)- However, for these high angular momentum modes to 
lead to annihilation, the gluinos must effectively tunnel through a large barrier. This means 
that direct annihilation is exponentially suppressed for high j. A significant annihilation 
cross section through larger angular momentum partial waves typically requires either that 
the object be of uniform density (i.e. like a macroscopic composite system) or that there are 
high angular momentum bound states that contribute to the cross section. 

One potential way to avoid the tunneling suppression at large angular momenta would 
be via the exchange of a high j QCD resonance. In this case, the resonance itself would 
carry the angular momentum, and could lead to immediate annihilation. While this might 
be plausible if the gluinos had a mass in the GeV range, it does not seem plausible for 
TeV mass gluinos- there are no sharp QCD resonances well above the GeV scale. If such a 
resonance were to exist, it would be extraordinarily broad, and would therefore not lead to 
rapid resonant annihilation. 

A second way to bleed off angular momentum would be by radiating. For the relevant 
large angular momentum states, we can use a semi-classical treatment where radiation is 
caused by accelerating a particle. This acceleration could be caused if a QCD string formed 
between the two gluinos. In this case, the radiation may be described via the Larmor formula 
where the power radiated is 

p ~ a2 ~ rr^2' ( 5 ) 

where a' is the tension of the exchanged QCD string. The total energy radiated can be 
estimated by multiplying this power by the crossing time t cross ~ R/v, where v represents 



3 




100 200 500 1000 2000 10000 
Gluino Mass (GeV) 

Figure 1: Gluino abundance per co- moving volume as a function of mass. Three curves are 
shown. In the first (solid), the annihilation cross section is assumed to be simply given by 
the perturbative cross section of Eqn. The other curves correspond to a cross section that 
saturates s-wave (dashed) and s-wave plus p-wave unitarity (dot-dashed). 



the relative velocity of the two hadrons. For T ~ Aqcd, we find 

E rad ^(^£R) 3/ \ QCD . (6) 
V m g / 

Thus, radiation from the gluinos is small, in fact, much less than the mass gap to the 
lightest possible state that could be radiated (the pion). This agrees with the intuition that 
heavy objects do not radiate. We should note that we could have applied a similar Larmor 
argument prior to the QCD phase transition. In this case, the relevant force is not due to 
a QCD string, but rather to a QCD Coulomb potential. In this case, the radiation will be 
further suppressed by perturbative powers of a s , again arguing against a large rate for the 
formation of bound states. Light QCD degrees of freedom do not carry the momentum or 
angular momentum of the system. Radiation from the cloud therefore is not able to reduce 
the relative angular momentum of the heavy gluinos, so they remain incapable of direct 
annihilation. 
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2.3 Relic Abundance Summary 

In summary, we estimate the total annihilation cross section as follows: for T > Aqcd, we 
simply utilize the perturbative cross section for gluinos jS]; for T < Aqcd, we allow for the 
possibility of an increased cross section. However, because we have argued that high angular 
momentum states do not significantly contribute to the annihilation of gluinos, we do not 
allow for an arbitrarily large cross section at this stage. We expect that the annihilation is 
conservatively given by a cross section that saturates s-wave unitarity. We take Aq C d = 200 
MeV as the point where the cross section changes. 

This cross section is thermally averaged ^U] and used in the numerical integration of the 
Boltzmann equation. We treat the QCD phase transition as second order although this is 
not expected to significantly affect the results [TT] . 

For completeness, in Fig. 2, we show the relic abundance for three cases. First, a solid 
curve that shows the relic abundance, assuming all annihilation is perturbative. Second, a 
dashed curve that incorporates the turn on of a cross section that saturates s-wave unitar- 
ity. Finally, a dot-dashed curve that has a cross section that saturates s-wave plus p-wave 
unitarity, which we view as an even more conservative assumption. 

3 Limits 

In this section we discuss the effects of the decays of the relic gluinos. Depending upon the 
lifetime of the gluino, these decays can disturb the predictions of BBN, distort the Cosmic 
Microwave Background Radiation (CMBR), or show up in the diffuse gamma ray back- 
ground. If the gluinos have a lifetime of order the age of the universe, they can potentially 
be observed in experiments searching for heavy exotic nuclei. We now discuss each of these 
constraints. 

First, we consider the era of BBN. These constraints are a function of the gluino de- 
structive power, £ = E vis Y. Here E vis is the energy of the gluino decay products which 
are deposited in the thermal bath, and Y = n/s represents the number density of gluinos 
per co-moving volume. To determine E^s, we assume that half of the energy of a decaying 
gluino goes into the neutralino, while the other half is distributed between hadrons. At the 
time of the gluino decay, the neutralino mean free path is large, and thus the energy in the 
neutralino is simply carried away; it is not "visible" . That is to say, we are taking the visible 
energy E vis = rag/2. If the lightest supersymmetric particle has a mass close to that of the 
gluino, the limits are weakened for a given abundance. 

Constraints arising from BBN on the lifetime of a hadronically decaying particle as a 
function of its destructive power were described in We take the most conservative D/H 
and 4 He limits. We are interested in the limits as a function of the gluino mass, but they 
only present results for three masses (100 GeV, 1 TeV, and 10 TeV). For a fixed lifetime, 
the dependence of £ on m g is roughly power-law, and so we can use these three points 
to interpolate to different values of the gluino mass. This procedure gives bounds on the 
gluino lifetime as a function of mass. We then convert lifetime bounds to bounds on the 
supersymmetry breaking scale [121 E], augmenting the decay rate by the two-body decay 
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rate when significant. The result is shown in Fig. 2. 

We now discuss in more detail the origin of the BBN bounds. Three BBN bounds are 
particularly restrictive: D/H, 6 Li/H, and 3 He/D. After BBN has finished, the universe is 
composed primarily of hydrogen and helium, with trace amounts of other elements. A gluino 
decay releases particles at energies much higher than the typical kinetic energies of nuclei at 
that time. At times later than about 100 seconds, mesons decay before they scatter, so the 
most destructive hadronic decay products are baryons. Protons produced by gluino decays 
scatter off of the background protons and 4 He nuclei. The elastic scattering increases the 
energy of background protons, making them more likely to react with other elements. The 
inelastic p + 4 He — ■> 3 He + D reaction increases D/H. If a D scatters off of a background 4 He 
nucleus, 6 Li can form. The epoch important for D formation is around 100 seconds - this 
sets the lower edge of the BBN curve in Fig. 2. 

A second process important in BBN is photo-dissociation ^H]- Photons with energies 
above 20 MeV can break 4 He nuclei into 3 He or Tritium (which later decays weakly to 3 He). 
This can cause the 3 He/D ratio to become too large. However, this process does not become 
important until the temperature drops to ~ 1 keV at ~ 10 6 seconds (see, e.g. [TB]). At 
temperatures before this, the 20 MeV photons lose energy efficiently by scattering off the 
background photons before they can break apart a nucleus. 

Gluinos that decay after recombination can give rise to photons that free-stream to us, 
and are visible in the diffuse gamma ray background JTj. Photons are produced when gluinos 
decay to pions which subsequently decay to gamma rays. Observations by EGRET [18 j and 
COMPTEL limit the flux of these gamma rays and thus the relic abundance of gluinos. A 
3-body decay including an invisible product, in our case the neutralino, was assumed in HT1. 
This is also the case at hand. The region excluded is shown in Fig. 2. These observations of 
the diffuse gamma ray background rule out a gluino with lifetime around the present age of 
the universe. 

The CMBR can be used to put further limits on the relic abundance of gluinos with 
lifetimes up to ~ 10 13 s [T^ I20j. Gluinos that decay during or after the epoch of thermal- 
ization of the CMBR can distort its spectrum and are thus limited by observations from 
COBE[2^. To derive these limits we make the conservative assumption that roughly 10% of 
the energy of the gluino gets transmitted to the CMBR when it decays, with the remainder 
carried off by neutrinos or the lightest supersymmetric particle. The limits derived from 
these constraints are shown in Fig. 2. 

If the gluino lives to the present day we would expect it to bind into nuclei, producing 
anomalously heavy isotopes. A combination of time of flight, mass spectrometry, and direct 
density measurements places severe limits on the abundance of terrestrial heavy elements 
today [22 l2"3j . This rules out a gluino with lifetime greater than the present age of the 
universe for the entire gluino mass range considered. These bounds are weakened if it is 
assumed that the heavy elements sink towards the center of the earth. However, for heavy 
hydrogen, the equilibrium time constant is > 10 8 years and so if the oceans undergo mixing 
the heavy hydrogen would be roughly uniformly distributed [22J. The bounds derived from 
these searches are shown in Fig. 2, labeled "Heavy Hydrogen." 

Gluino properties are already restricted via direct searches at colliders |24l 12*5] . The 
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Figure 2: Limits on the supersymmetry breaking scale, m>s, as a function of the gluino 
mass, rrig. The bounds are derived assuming a perturbative cross section in Eqn. |3] for 
temperatures greater than the QCD phase transition, T > 200 MeV. For T < 200 MeV, we 
assume that the annihilation cross section saturates s-wave unitarity. Also shown (dashed) 
are the limits in the case where the annihilation cross section saturated s-wave plus p-w&ve 
unitarity. The shaded regions are excluded. The lower edge of the BBN curve arises from 
the requirement that the D/H ratio remained undisturbed, and corresponds to a lifetime of 
approximately 100 seconds. 



weakest limits come in the case where the produced i?-hadrons are neutral. Then, they 
escape the detector, carrying away a substantial fraction of the event energy. The gluinos 
are then observed in mono-jet events, triggered by the presence of an additional radiated jet. 
Current limits on the gluino from Tevatron Run I are m~ g > 170 GeV [25] • This bound should 
increase to 210 GeV if Run II sees nothing, and to 1100 GeV at the Large Hadron Collider 
(LHC). These bounds are independent of the supersymmetry breaking scale. If charged 
i?-hadrons are produced that do not immediately decay, this bound improves. Searching 
for gluinos through anomalously slow tracks in the tracking chambers of both ATLAS and 
Tevatron were studied in [2SJI2ZI- These provide search reaches comparable to the monojet 
signal and will be useful as an additional discovery channel. 

Finally there is the possibility of seeing gluinos in cosmic rays and was studied by [28j. 
If a gluino were seen, then this would set and lower limit on the susy breaking scale. While 
most of the available lifetime and mass ranges are ruled out by the considerations in this 
paper, there is a window available at low gluino masses for gluinos to be seen in cosmic rays. 



7 



4 Conclusion 



In this note we have explored the limits on the split supersymmetry parameter space coming 
from cosmological constraints on a long lived gluino. If the gluino were to annihilate with 
a geometric cross section after the QCD phase transition, then it would evade most cosmo- 
logical constraints. However, we have argued that this is not the case, the cross section is 
more likely set by the de Broglie wavelength of the gluino. Then, the cross section is likely 
bounded by saturating s-wave unitarity. Using this cross section, it is possible to extract 
useful limits from the CMBR, BBN, and the diffuse gamma ray background. For gluinos 
heavier than 300 GeV the earliest cosmological constraints are from BBN, specifically the 
D/H ratio. This sets an upper bound to the lifetime of roughly 100 seconds and an upper 
bound on iris of 10 9 GeV. For gluinos lighter than 300 GeV, the earliest constraints arise 
from non-observation of diffuse gamma rays and sets an upper limit to the lifetime of 10 6 
years and corresponds to ms < 10 11 GeV. If collider searches at the LHC find a gluino 
heavier than 600 GeV, then the lifetime will be less than 100 seconds. This has implications 
for experiments that hope to trap the gluino and directly measure its lifetime. 
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